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Abstract 

Cued and un-cued cell polarization is a fundamental mechanism in cell biology. As an alternative 
to the classical Turing bifurcation, it has been proposed that the cell polarity might onset by means 
of the well-known phenomenon of wave-pinning (Gamba et al, PNAS, 2005). A particularly simple 
and elegant model of wave-pinning has been proposed by Edelstein-Keshet and coworkers (Biop. 
J., 2008). 

However, biomolecular networks do communicate with other networks as well as with the external 
world. As such, their dynamics has to be considered as perturbed by extrinsic noises. These noises 
may have both a spatial and a temporal correlation, but any case they must be bounded to preserve 
the biological meaningfulness of the perturbed parameters. 

Here we numerically show that the inclusion of external spatio-temporal bounded perturbations 
may sometime destroy the polarized state. The polarization loss depends on both the extent of 
temporal and spatial correlations, and on the kind of adopted noise. 

Namely, independently of the specific model of noise, an increase of the spatial correlation induces 
an increase of the probability of polarization. However, if the noise is spatially homogeneous then 
the polarization is lost in the majority of cases. 

On the contrary, an increase of the temporal autocorrelation of the noise induces an effect that 
depends on the noise model. 
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I. INTRODUCTION 



The formation of two distinct spatial domains able to distinct in a single cell two sep- 
arated parts, let us say a "tail" and a "head", is called cellular polarization [T]. This 
phenomenon is at the basis of two fundamental bioprocesses: the asymmetric cell division 
and the chemotactic motion of some kinds of cells (for example neutrophils) [I]. 

The first is a patterning that has to regularly repeat until the replicating potential of a 
proliferating cell is exhausted (theoretically never for stem cells), thus it is thought to be 
essentially un-cued However, there are models, such as that by Ortoleva and Ross [3], 
that postulate a role of random cues also for the onset of asymmetric division. Finally, also 
in absence of mitosis, a cell can experience symmetry-breaking and start moving in random 
directions 

On the contrary, the patterning during chemotactic motion induced by an external cue: 
the gradient of the chemoattractant to be followed (or of the chemorepellor to be avoided) 

urn- 

In both cases it is necessary that a pattern "head" -"tail" must be formed. From the bio- 
physical point of view, both cases are classified as pattern onset in non-equilibrium systems. 
Patterns formation in biosystems is a fundamental topic of computational biology, and it 
is largely influential in experimental biology. This area of investigation has been started 
by Alan Turing's very well-known paper on morphogenesis |7j, where he modeled the onset 
of a pattern in two kinds of multi-cellular structures (a ring and a sphere) as a symmetry- 
breaking bifurcation driven by a strong difference in the diffusion of two morphogens. This 
mechanism is called the Turing's bifurcation. In early seventies the Turing mechanism has 
been biologically substantiated by Gierer and Meinhardt [8HTO] , who introduced two possible 
patterning mechanisms based on Turing bifurcation. The first is the Activation-Inhibitor 
model, where the pattern is induced by the reaction-diffusion interplay between a short range 
self-activator (which is membrane-bound or, any case, has a very low diffusion coefficient) 
and a long-range inhibitor (which has a large diffusion coefficient). The second model is 
the Activation-Depletion model, where two proteins interact: the first is again a short-range 
self-activator, but the second is no more an inhibitor. On the contrary, the second protein 
is, instead, depleted by the activator. 

These pioneering studies generated a wide amount of literature, which is increasingly 



2 



linked with experimental data [21 [TlHT^ . 

In 2005, Gamba and colleagues [T3HTS] proposed a stochastic model of chemotaxis-induced 
spatial symmetry breaking in a single cell. The central core of this model was the interplay 
between two membrane-bound molecules and two cytosolic molecules. The presence of some 
feedbacks (with related bistability) and the difference of diffusion coefficients induced a 
mechanism of phase separation that is different from the ones based on Turing's bifurcation. 
Quite interestingly, in the supplementary materials of [15] Gamba et al. proposed a simplified 
and mean-field version of their model that exhibits, in one spatial dimension, a traveling 
wave whose velocity decreases until the front stops. 

In 2008, Edelstein-Keshet and colleagues proposed in [T9H22] a similar but simpler 
biomechanism leading to the onset of single-cell polarization due to the 'freezing' of an 
initially traveling wave of large concentration of a membrane-bound protein. Indeed, this 
elegant and minimalist mechanism is based on the inter-conversion of a membrane-bound 
active protein 'A' in its cytosolic inactive form 'B', where 'A' positively feedbacks on its 
activation. Also this process induces the initial onset of a wave-front that stops. This phe- 
nomenon is called in physics 'wave-pinning'. The wave- pinning- induced polarization will be 
here called WPP. 

The WPP process is biologically back-grounded, since it well fits to some key pro- 
teins involved in cell polarization [T9|, such as Rho-GTPases, which switch between active 
membrane-bound and inactive cytosolic forms. 

Some similarities exists between the WPP and the Activation-Depletion induced polar- 
ization (ADP) in that: i) in WPP the two involved molecular forms have very different 
diffusion coefficients; and ii) the positive feedback of 'A' on its activation that depletes 'B' 
may remember the depletion mechanism. However, in WPP the patterning mechanism is 
totally different, since in it no Turing bifurcation is observed, and, conversely, in the Turing- 
based models no waves are observed. Moreover, one of the 'essential ingredients' [IH] of 
WPP model is the presence of bistability (see also [IHl [IT]), making it a genuinely non- 
linear model, whereas Turing bifurcation stems from the linearization of reaction-diffusion 
equations. Finally, in WPP the total mass of the two forms 'A' and 'B' is conserved. 

Summarizing, from a theoretical point of view, although they are two quite different 
mechanisms, both ADP and WPP have equivalent effects of inducing a pattern. However, 
passing from the general principles to the concrete numerical simulations in case of realistic 
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values of the parameters, in [19] it has been shown that the WPP has (at least for the 
examined range of parametric values) an important advantage on ADP: it is extremely 
rapid. Indeed, in WPP the pre-patterning transient has a characteristic time of more or less 
10 seconds, whereas ADP would require more than 10 minutes [19]. The latter time is too 
long and would make ADP unsuitable, especially for chemotaxis-driven polarization. 

This temporal difference in the duration of respective transients makes WPP a very 
attractive mechanism to be experimentally validated. 

Moreover, from the point of view of theoretical physics, the WPP model is one of few 
examples (see also [23]) of robust wave-pinning. Indeed, standard wave-pinning condition is 
fulfilled for isolated points of a key parameter, and also for very small parametric changes 
the wave restarts traveling. 

Very recently, it has been published a paper [20] on the robustness of the WPP model 
with respect to the intrinsic stochasticity, i.e. by considering the fact that the molecular 
distribution is discrete and not continuous. The main result of that investigation is that for 
very small number of molecules the wave collapses and the WPP is destroyed. 

Here we want to study another an equally important problem: the robustness of WPP 
with respect to the unavoidable presence of extrinsic noise. 

As we mentioned above, one of the basis of WPP is the presence of multistability. The 
interplay between extrinsic noise and multistability is of fundamental relevance in chemistry 
and systems biology [21H2S] , as shown since the pioneering work by Kramers [2Z] . In a non- 
linear system, indeed, the presence of noise may cause the transition from a state to another 
[28] (and other important effects [26]). This in the WPP case might mean that the noise 
could induce the switch from a polarized to a non-polar molecular distribution. Moreover, 
the presence of noise may deeply affect the dynamics of traveling waves, as reviewed in [29] . 

Moreover, the above-mentioned rapid duration of the transient leading to the full po- 
larization suggests that the WPP might not be robust to realistic colored extrinsic noises 
representing the interplay of the WPP mini-network with other biomolecular networks that 
have the same or even larger characteristic times. 

We investigate here this point by means of numerical simulations of the WPP model in 
the case that its parameters are perturbed by bounded colored spatio-temporal noises. 
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II. UNBOUNDED AND BOUNDED NOISES 



Let us consider the well-known zero-dimensional Ornstein-Uhlenbeck stochastic differen- 
tial equation: 

c{t) = —m + —v{t), (1) 

Tc Tc 

where Tc > 0, \/2D is the noise strength and 77 (t) is a Gaussian white noise of unitary 
mtensity{T](t)r]{ti)) = 5{t — ti). It is well-known that solution of equation ([T| is a Gaussian 
colored stochastic process with autocorrelation: (^(t)^(ti)) oc exp{~\t — ti\/Tc). In [30] 
equation ([T]) was generalized in a spatially extended setting by including in it the most 
known and simple spatial coupling, the Laplace operator, yielding the following partial 
differential Langevin equation: 

1 \/2D 

dt^x, t) = a^oise^^^x, t) - t) + r^{x, t), (2) 

Tc Tc 

where Dnoise > is the spatial correlation strength [30] of ^(x, t). 

Two strategies have mainly been adopted so far for defining zero-dimensional bounded 
noises: i) nonlinear filtering unbounded noises; ii) applying bounded functions to unbounded 
noises. 

The first approach was employed in [3T1433] , where the following family of bounded noises 
was introduced: 

«W = -i«W + /J^.W. (3) 

where rj{t) is a Gaussian white noise with {rj{t)rj{ti)) = 6(t — ti). The bounded nature of the 
noise described in ^ easily follows from the fact that ^ = +B {—B) implies ^ < (> 0). 
The process S,{t) has zero mean and the same autocorrelation of the OU process [31l 132] , 
and its stationary probability density is given by: Pcl(0 = ^ (-^^ " + - For ^ > the 
distribution is unimodal and centered in 0, while for —1 < 2; < it is bimodal, having a 
"horned" distribution with two vertical asymptotes at ^ — )■ ±5. 

In order to define spatio-temporal bounded noises based on the Cai-Lin noises, in [31] 
we adopted an approach analogous to the one employed in [301 EH] and in equation ^ to 
extend the OU process. Thus, the spatio-temporal Cai-Lin noise can be defined as follows: 



dti{x, t) = D^oiseV'^ix, t) - t) + J , t). (4) 

Tc V Tc[l + Z) 



The sine- Wiener noise [36j is obtained by applying the bounded function h{u) 



B sm[-\/2/TcU) to a random walk W{t) defined as W = //(t), where r]{t) is a white noise of 
unitary intensity, yielding: 

m = Bsin{J^W{t)^. (5) 

The stationary probability density of C,{t) is given by Psw{0 = V-^^ ~ C^); thus: 
Peq{±B) = +00. 

In [37], as a natural spatial extension of the sine- Wiener noise, we defined the following 
spatio-temporal noise: 

C(x,t) =5sin(27r^(a;,t)), (6) 

where $,{x,t) is defined by equation (|2|. 

Detailed studies of the properties of the above-defined spatio-temporal noises have been 
performed in [3ll [37] . 

III. BACKGROUND ON THE WPP MODEL 

The WPP model describes the interplay of two different forms (active, 'A', and inactive, 
'B') of a biomolecule, of which the active one is membrane-bound and the other is located 
in the cytosol and has a very large diffusion coefficient. Denoting with a{x, t) and b{x, t) the 
concentrations of, respectively, A and B, the WPP model reads : 

dta = D,V'a + f{a,b;p) (7) 
dtb = D,V%-f{a,b;p), (8) 

where Da and are the diffusion coefficients of A and B, and p is a vector of parameters 
(see below). Since A is membrane-bound and B is free in the cytosol, the following constraint 
must hold: 

Dt > Da. (9) 

A central hypothesis in [19j is that A positively feedbacks on its activation. Thus, it is 
assumed that /(a, 6;p) is of the form g{a;p)b with dag{a) > 0. Namely, in [19j it is employed 
the following functional form: 

/(a, b;p) = b (ka + j^^^^) - 6a. (10) 



As a consequence, the vector of parameters reads 

p = {ka,'y,K,u,6). (11) 
The boundary conditions are assumed of the no-flux type: 

Va(t, 0) = Va(t, L) = Vb{t, 0) = V6(t, L) = 0, (12) 
which, of course, imply the mass conservation: 

L 

{a{t,x) + b{t,x))dx = Q. (13) 

Note that the fundamental assumption ^ Da ensures the correct working, and the 
robustness, of the wave-pinning mechanism. 

In [in], Mori and co-workers showed that not only the above system is able to generate 
a polarization in response to initial transient cues and in absence of cues, but also that the 
onset of this pattern is faster than in the Turing mechanism. Indeed, they showed that the 
time to polarize is of the order of 1 to 10 seconds, and the effective time to 'complete the 
polarization' [19] is of about 30 seconds. 



IV. INCLUDING EXTRINSIC NOISE IN THE WPP MODEL 

We phenomenologically take into the account the interplay with extrinsic noise by as- 
suming, in the most general case, that: 

dta = DaV^a + f{a, b; p{x, t)) + ks{x, t)b (14) 
dtb = D,V% - f{a, b; p{x, t)) - ks{x, t)b, (15) 

where ks{x,t)b is the initial transient cue, which is null after a short time, and: 

Pi{x,t)=p^il + Ux,t))>0, (16) 

where ^i{x, t) is a bounded noise of the Cai-Lin or sine- Wiener type. We employ two different 
kinds of noises in line with the recent literature on bounded noises [Ml [371439] . which in 
other context showed that the statistical characteristics of a system perturbed by a bounded 
stochastic process depends not only on the bound of the noise, but also on its finer structure. 
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Unfortunately, no analytical tools are currently available to investigate the effects of 
bounded spatio-temporal perturbations, so that we shall resort to numerical simulations. 

We consider two kind of numerical experiments corresponding to, respectively, cued and 
un-cued polarization of a cell: 

i) external graded transient cue, modeled by: 

ks{x,t) = s{t){l-x/L) (17) 

S te[0,ti\ 
'5(1-1^) telhM- (18) 

T>t2 

a) random initial conditions with no external transient cue, i.e.: 

a{x, 0) = R r/(x)a_, b{x, 0) = 2.0, (19) 

where a_ = 0.2683312 corresponds to the lower homogeneous steady state of WPP model 
and ri{x) is a spatial noise uniformly distributed in [0, 1]. 

In both cases Mori et al. [12] showed that the transitory lasts less than 200 seconds, 
at which time the transitory is well established. However, the presence of spatio-temporal 
correlations might induce slower transients, so that we simulated all systems up to time 
T = 600 s. A first statistics we measured for each realization of our stochastic process is 
the following: 

W = i^'/^-J^f . (20) 

Jq a{x,T)dx 

which is the average 'position' x weighted by the normalized version of a{x, T). This statistics 
measures where the polarization is oriented. Of course, in the case of cued polarization this 
measure is mainly relevant to assess whether the polarization is in line with the cue, or if 
it is not present in the final simulation time. The second and more important statistics 
should then measure some variance of the density. One might be tempted to measure also 
this statistics in the final simulation time. However, during our numerical simulations we 
observed cases where the stochastically perturbed front temporarily became flat to then 
recover a polarized state, and so on. As a consequence, the analysis of the variability of the 
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profile of a must also consider the past history of the system. Namely, starting from the 
instantaneous amplitude of the distribution: 



St = (Max^e[ox]a(a;,t)) - (Max^e[o,L]a(x, t)) 



(21) 



we defined the following statistics: 



A 



t = niintg[2oo,T]'^t 



(22) 



Note that in case of noisy initial conditions, we expect an equiprobable distribution of the 
polarization, and also cases where both At is 'large' and (x) ~ 0.5, which indicate the 
presence of a 'central hump' in the distribution of a. 

As far as the characteristics of the employed noises are concerned, we shall assume that 
their temporal autocorrelation has a time-scale of 10 seconds. Also far larger values would 
be acceptable, since there is a number of biomolecular processes that are characterized by 
long time-scales. However, in those cases the loss of polarization would be almost sure. As 
far as the spatial scales of the noises, we considered three cases: i) spatially white noises; 
ii) spatially uniform noises; and in) finite spatial correlation. 

V. NUMERICAL SIMULATIONS 

In line with [19], in all our simulations we set the following values for the parameters 
of the WPP model: Da = 0.1 fim^s-\ Db = 10 fim^s-\ u = 1, 6 = 1 s-\ 7 = 1 s-\ 
ka = 0.067 s-\ K = 1 and cell size L = 10 fxm [HI HQHIS]. 

In all our simulations, the stochastic predictor-corrector algorithm [43j for Ito partial 
stochastic differential equations was adopted to simulate the noise. Moreover, in case of the 
Cai-Lin noise the transformation described in [34] was also adopted. The WPP model was 
then simulated by means of the second-order Runge-Kutta algorithm. 

Finally, a remark: in next sections, for the sake of the simplicity, instead of writing 
phrases like "this simulation of the model suggests that the cell does something" we shall 
concisely write (with an abuse of meaning) " the cell does something" . 
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A. External Cues 



In the first round of simulations we stimulated the system with the graded external 



cue, described in equations (17)-(18), with S = 0.07m s^^, ti = 20s and t2 = 25s. The 
initial conditions were: {a{x, 0),b{x, 0)) = (a_, Bq) ^ (0.2683, 20), corresponding to the lower 
homogeneous steady state of the unperturbed system [19] . In absence of noise we obtained, 
as in [19j, that the cell polarizes with time-scales of the order of 50s. 

By introducing an extrinsic noise in one of the relevant parameters of function /(a, b) via 



formula (16), we obtained a complex pattern of responses, in dependence of the spatial and 
temporal parameters of the noise, as well as of its type. 

In some cases, which we are going to illustrate, our stochastic model predicts that the 
cell polarization is preserved, although, of course, the wave-front experiences some random 
fluctuations; in other cases the cell depolarizes. In these cases, the profile of the concentration 
of A is fiat and oscillating. 

We never observed an inversion of the polarization, i.e. the polarized state follows in all 
cases the external cue. 

Typical configurations are shown in figure [l| where both a scatterplot (At, (x)) and three 
specific realizations are shown. 

We now illustrate the statistical response of the model to bounded perturbations that are 
spatially white, i.e. Dnoise = 0. 

Initially we shall consider perturbations in the parameter 6. 

In such a case, we observed that the response to the noise strongly depends on the type 
of the perturbation. Indeed, setting Tc = 10 s and B = 0.2, both Cai-Lin and sine- Wiener 
noises cause in many cases depolarization of the cell. This is illustrated by the scatterplots 
shown in figure [2| and in the distribution of shown in the upper panels of figure |4j 

By increasing the temporal autocorrelation up to 50 s we observed two dichotomic 
behaviors. If the perturbation is a Cai-Lin noise then the increase of causes a larger 
number of depolarizations (see figures [3^ and [4]:). On the contrary, if the bounded noise 
is of sine- Wiener type then one observes a larger probability that the cell maintains the 
polarization induced by the external deterministic cue (see figures [3|d andjiji). 

As far as perturbations in u and 7 are concerned, fluctuations in 7 induce effects com- 
parable to the ones caused by noises affecting S, whereas perturbations of u depolarize to a 
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lesser extent. Note that the perturbation of oj may summarize the effect of extrinsic noise 
in the feedback mechanism. 

We now examine the effects of spatially correlated noises. In figures [6]j7] we show that the 
presence of non-null spatial correlation, i.e. Dnoise > 0, increases not only the probability of 
maintaining the cell polarization, but also its intensity (i.e. the magnitude of spatial gradient 
of a). Indeed, in such cases simulations suggest not only a decrease of the probability of 
observing small or null A^, but also a significant increase of the probability that is 'large'. 
For example: if Dnoise = 0.01 the probability -P(O) is quite large (in 30% of cases the cell 
is not polarized), whereas for Dnoise > 0.1 -P(O) is small and it is a decreasing function of 
D 




(a) 




(b) 





(c) 



(d) 



FIG. 1: Bounded stochastic perturbation of parameter 5. Type of noise: spatially white Cai-Lin 



noise with temporal correlation Tc = 10 s. Panel (a): scatterplot (A^, (x)) for the density a{x,t). 



For large and intermediate values of the system is polarized, see panels (c) and |(d)[ and the 
fluctuations of the polarization front are more pronounced for decreasing At- For A^ ~ the system 
is not polarized, and is characterized by a (spatially constant) temporal oscillation dynamic, see 
panel |(b)[ 
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FIG. 2: Spatially white but temporally colored perturbations of parameter 5: scatterplot (A^, (x)). 
Panel (a) Cai-Lin noise; panel |(b)[ sine- Wiener noise. In both cases: Tc = 10 s. Number of 
simulated systems: 50. 
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(b) 



FIG. 3: Spatially white but temporally colored perturbations of parameter 5: scatterplot (A^, (x)). 
Panel [(a)| Cai-Lin noise; panel [(b)| sine- Wiener noise. In both cases: Tc = 50 s, B = 0.2. Number 
of simulated systems: 50. 

We also simulated the response of the WPP model to a noise that is spatially uniform 
but temporally varying. Namely: i) instead of the spatio-temporal Cai-Lin noise, we used a 
noise P{x,t) = ^{t), where C,{t) is the temporal Cai-Lin noise; ii) instead of the sine- Wiener 
noise we employed a noise f3{x,t) = Bsin{2n^(t)) where ^(t) si the OU noise. Since in this 
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FIG. 4: Spatially white but temporally colored perturbations of parameter 6: heuristic histograms 
of the distributions P{At). Upper panels: Tc = 10 s; lower panels: Tc = 50 s. Panels [(a)] and 
(c) response to Cai-Lin noise; panels |(b)| and |(d)[ response to sine- Wiener noise. In all panels: 
B = 0.2. The response to increasing Tc is different for the two kinds of noises: the depolarization 
probability increases for Cai-Lin noise, decreases for sine- Wiener noise. 



case the spatial correlation length is infinite, we expected to observe that all the simulated 
cells could maintain their polarization. Quite interestingly, we observed that this spatially 
uniform noise induces (see fig. [7]F) in the model a behavior that is comparable to the one 
induced by spatially white noises (compare with previous fig. [7k). 
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FIG. 5: Spatially white but temporally colored perturbations affecting uj (panel (a) ) and 7 (panel 



|(b)[ ): heuristic histograms of the distributions P{At). In both cases: Cai-Lin noise with Tc = 10 s 
and B = 0.2. 
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FIG. 6: Spatially and temporally colored noise: scatterplot (A^, (x)) for various values of Dnoise- 
Here the perturbed parameter is 7. Type of noise: Cai-Lin. Parameters Tc = W s, B = 0.2. 
Number of points for each series is 200. 



14 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
A, 



(a) 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 ^ 
0.1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
A, 



(c) 



0.9 
0.8 ^ 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
A, 



(e) 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




D„„i5P=0.01 




(b) 



0.9 




0.8 




0.7 




0.6 




0.5 




0.4 




0.3 




0.2 




0.1 







n 






0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



(d) 



0.9 
0.8 
0.7 

0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




uniform sp. corr. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
A, 



(f) 



FIG. 7: Histograms of P{^t) in case of various degrees of spatial correlation of a Cai-Lin per 
turbation applied to parameter 6. Panel (a) Dnoise = 0; panel |(b)[ D 
Dnoise = 0.1; panel [(d)] Dnoise = 1; panel [(e)J Dnois, 
Other parameters: Tc = 10 s, B = 0.2. 



; = 0.01; panel (c) 
10; panel |(f)[ spatially uniform noise 
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B. Initial random state 



Here we illustrate the results of simulations performed without adding the external de- 
terministic cue, but considering an initial random distribution of A , following the equation 



([19j) with R = 3. 

In absence of noise, we observed three possible kinds of steady states: polarization (left 
or right, in this case it does not matter) or a central 'hump'. 

As far as spatially white noises are concerned, we observed a behavior depending on the 
perturbed parameter, but not on the type of noise. Indeed: 

• in case of perturbation of S (see figure |8| where the noise is of sine- Wiener type): for 
B = 0.2 solutions with central hump are not observed, and in some cases cell is not 
polarized. Absence of polarization is always observed if B = 0.35; 



• in case of perturbation of a; (see figure |9]) or of 7 (not shown): for B = 0.2 cells seldom 
is not polarized and the number of solutions with central 'hump' is roughly similar to 
that observed if B = 0.05. For B = 0.35 in sometime cases the cells remain polarized 
(and no cells has the 'hump' pattern). 

Finally, also here the increase of spatial correlation restores the polarization. 



VI. CONCLUDING REMARKS 



Our main aim in this work was to study the robustness to extrinsic stochasticity of an 
interesting mechanism explaining both cued and un-cued cellular polarization: the wave- 
pinning [19] . 

We modeled the above-mentioned extrinsic perturbations by means of spatio-temporal 
bounded noises that model perturbations in a more realistic way than Gaussian noises, which 
are unbounded. In particular, we adopted the following two noises: the Cai-Lin [31] and the 
sine- Wiener [37] noises, which we recently proposed. 

The robustness of the WPP model to extrinsic noise was, in particular, investigated by 
means of suitable statistics summarizing the polarization probability, the fluctuations and 
the intensity of the polarization. 

An interesting feature of the deterministic WPP model is that it is deterministically 
robust, in the sense that most models exhibiting the wave-pinning phenomenon are such 
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FIG. 8: Random initial conditions and spatially white sine- Wiener perturbation of 5. Scatterplot 
(At, (x)) for three values of B. Other parameters: Tc = 10 s. 

that the waves 'freezes' for a specific isolated value of a given parameter [T9l 123] . In the 
WPP model, instead, the system self regulates and the wave freezes for an entire range of 
the chosen bifurcation parameter |19j. 

Here we showed that under extrinsic stochastic perturbations, the maintenance or loss of 
polarization strongly depends on both spatial and temporal colors, as well as on the specific 
kind of noise. 

Indeed, in case of Cai-Lin noise the increase of temporal correlation induces a decrease 
of polarization probability, whereas in case of sine- Wiener noise one observes the opposite 
phenomenon. 

Moreover, in both cases passing from white to spatially correlated noise induces a loss of 
polarization. However, in the simulations where we employed a spatially constant noise we 
again obtained a large probability of polarization. 

In other words, the WPP model exhibits a sort of contrast of colors (the spatial and the 
temporal) and of noise types. 

Concerning the possible use of the WPP model to describe the spontaneous cellular polar- 
ization emerging from random initial conditions, beside the polarized and global oscillating 
states, we observed (in line with (TH]) the emergence of "non-polar" patterns that are char- 
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<x> 

FIG. 9: Random initial conditions and spatially white Cai-Lin perturbation of u. Scatterplot 
{At, {x)) and corresponding histograms of the distributions P{At) and P({x)). Other parameters: 
Tc = 10 s. 

acterized by a hump in the density of A, located in the center of the cell. Moreover, the 
effect of the amplitude of the noise (and of the spatial coupling) is 'non- monotone'. Indeed, 
it is possible to find an intermediate amplitude of the noise (or spatial coupling strength) 
by which both the probabilities of unpolarized oscillating states and of 'humped' states are 
low. 

Summarizing the above observations, we may say that the answer to the main substantive 
question 'Is the WPP model robust?' is: 'it depends on the context where the cell is 
embedded'. In other words, being equal the amplitude of the external perturbations, the 
robustness of the wave-pinning-based mechanism of cellular polarity strongly depends on 
the kind of extrinsic stochasticity that affects the WPP network. Note also that the ability 
of this mechanism to describe a spontaneous cellular polarization is questionable (as also 
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stressed in [21]) because of the presence of humped solutions. Although in presence of 
extrinsic noise the hump may disappear, however in those cases the onset of noise-induced 
spatially homogenous states is observed. 

As far as the onset of wave-pinning-induced polarization in case of low number of 
molecules, studied in [20], we remind the reader that for moderate number of molecules 
a stochastic reaction system may be approximated by the corresponding deterministic sys- 
tem plus suitable multiplicative external white Gaussian noises. This is a well-known result 
in statistical and chemical physics [HI |13], also known in Systems Biology as 'chemical 
Langevin equation' method [IH]- Thus here we might be tempted to read the result of paper 
[20] as fragility of the WPP mechanism to white noise perturbations. However, we stress 
here that for computational reasons the simulations of [20] has been performed by assuming 
that the number of molecules of the species B is spatially homogeneous. Moreover we recall 
that it is the equation ([s]) that guarantees the robustness of WPP model. This might imply 
that the stochastic version of WPP model could be more robust than it has been reported 
in the numerical simulations of [20] . 

Thus, our aim is to explore in a future work the full discrete stochastic version of equa- 
tions ((7])-([8]), also including extrinsic noises. Indeed, the interplay of intrinsic and bounded 
extrinsic noise may be of interest in Systems Biology, as stressed in [26] . 
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